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Approximation by Spline Curves: 
towards an Application to Cognitive Neuroscience 

Maria-Laura Tonente, Stefano Anzellotti, Chiaia Finocchiaro and Claudio Fontanari 


Abstract —We present a procedure to approximate a plane 
contour by piecewise polynomial functions, depending on various 
parameters, such as degree, number of local patches, selection 
of knots. This procedure aims to be adopted to study how 
information about shape is represented. 

Index Terms —Shape recognition, perceptual discrimination, 
spline curves, B-splines, least squares approximation. 


I. Introduction 

The shape of an object is a fundamental source of infor¬ 
mation for recognition. While color and texture information 
are also important (see ID), line drawings, which discard 
texture and color, are usually sufficient to recognize an object. 
Furthermore, infants generalize object labels on the basis 
of shape similarity rather than color similarity (see Gl), 
suggesting that shape plays a key role in categorization. In 
order to understand how we recognize and categorize objects, 
therefore, it is crucial to study how information about shape 
is represented. 

Since the space of shapes is infinite dimensional, the brain 
likely approximates it using a lower number of dimensions in 
order to make the problem more tractable. For comparison, 
it can be useful to consider the example of color. Like 
the space of shapes, the space of light spectra is infinite 
dimensional. However, with three dimensions it is possible to 
model accurately the space of perceptually discernible colors 
(see 0 ). In the same vein, we can look for finite dimensional 
spaces that encompass the perceptually discernible shapes. 

Recent work (see J4)) found that a model with approxi¬ 
mately 47 dimensions explains accurately a large amount of 
the variance in the object recognition errors and similarity 
judgments of human participants. However, models of this 
type leave open the question of whether the dimensionality of 
shape space is the same in the neighborhoods of all shapes. In 
addition, the most relevant dimensions for recognition might 
differ in the neighborhood of different shapes. For example, 
the relevant dimensions for discriminating between individual 
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faces and for discriminating between buildings might not be 
the same, and thus different low-dimensional approximations 
of the space of shapes could be used in the two neigh¬ 
borhoods. This possibility is compatible with the differential 
involvement of different brain regions within temporal cortex 
in the processing of different object categories (see 0 ). Our 
interdisciplinary research project aims to introduce techniques 
to individuate lower dimensional subspaces that locally ap¬ 
proximate the space of perceptually discernible shapes. 

A couple of crucial remarks are in order here. First, our 
approach aims to model shape discriminability rather than 
object recognition accuracy. In other words, it is a model of 
perception rather than a model of perceptual categorization 
(in this respect, it is quite similar to the case of color space). 
The choice of investigating discriminability is based on the 
idea that a common representation of shape underlies both 
shape discrimination and judgements of perceptual similarity 
that are at the basis of perceptual categorization. This idea is 
supported by the evidence currently available (see 0 ). Next, 
our goal is to represent shapes locally rather than constructing 
a single lower dimensional space underlying the representation 
of all shapes. This choice is a consequence of the aim to 
model perceptual discriminability: different local spaces seem 
necessary to account for the large variety of shapes we can 
perceptually discriminate. 

For modeling shape space in the neighborhood of a fixed 
shape spline curves will be used. Spline curves can ap¬ 
proximate contours closely, they enjoy elegant mathematical 
properties, and their computational complexity is relatively 
low. For all these reasons, splines are widely considered to 
be an ideal tool to approximate signals (see for instance ©)■ 
In particular, the use of B-splines for digital signal processing 
has a long history (at least since (TJ). In this paper we employ 
spline curves to address the following approximation problem: 
given a plane contour C construct a (parametric) spline curve C 
of a fixed degree d well approximating C. The procedure 
(described in Section m and based on Section E is divided 
into three parts: use of strengthened edge detection techniques 
to extract the points of the contour C; computation of the spline 
curve C, performed by solving a Least Squares Approximation 
problem (see Problem |II.5| and subsequent discussion); plot 
of the computed spline curve (and, eventually, comparison 
with C). The spline approximation part (second part) is the 
core of the presented procedure and, as it is well known to 
the experts (see for instance 0), an important issue for it is 
the choice of knots. At the end of Section [II] we recall three of 
the most used parametrization methods (uniform, chord length 
and centripetal). The chord length parametrization method is 
then chosen to show a concrete example (see figures [TTa)| 1(f) )>. 
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II. Basics of Spline Approximation 

The main reference for this section is 0- 
Bezier curves and, after their introduction in 1946 (see 
ESI). spline curves have been widely used to construct smooth 
curves from a given set of points in an efficient and numeri¬ 
cally stable way. Their geometrical construction is essentially 
based on recursive convex combinations of curves of smaller 
degrees. Namely: 


Definition II. 1. Let n,d be positive integers, with n > d + 1. 
Let pi,... ,p n £ R 2 be n control points, and let r = 
be the knot vector, which is assumed to be a nondecreasing 
sequence of n + d — 1 real numbers, that is, < ... < t ra +d- 
The functions p ij(t) are recursively defined in the following 
way: we set Pi j o(^) = Pi, for i = 1,..., n and 


P = tl+d 3+1 _ * Pi-i,j-i(*) + 

“i+d-j-j -1 £i 


t-ti 


ti+d—j+1 t>'i 


■Pi,j-i(*) 


for j = 1 ,d and i = j + 1,...,n, where possible division 
by zero is resolved by the convention that a division by zero 
is zero. The (parametric) spline curve f(t) of degree d with 
control points (pi)" = i and knots t = is defined as 


n 

fit) = PiAt) B i,o(t) 

i=d -\-1 


where Biflit) are the piecewise constant functions defined by 




1 if ti < t < t i+ 1 

0 otherwise 


Due to their special construction, spline curves enjoy many 
important geometrical properties, such as local dependence on 
the control points, being contained in the control polygon, and 
so on. 

A crucial fact is that every spline curve f(t) can be written 
as a linear combination of special univariate polynomials 
(the B-splines, see Definition II.2| ) with the control points as 
coefficients (see Theorem |II.3 >. 


Definition II.2. Let m, d be nonnegative integers, with m > 
d+ 2, and let r = (L)™, be a nondecreasing sequence of real 
numbers. Let i = 1,..., m+d+ 1; the ith B-spline of degree d 
with knots r(ti)2 =1 , denoted by Biflit), is recursively defined 
in the following way: for i = 1,..., m — 1, we set 


Biflit) 


1 if U < t < t i+ 1 

0 otherwise 


Then, for each j = 1,..., d and i = 1,..., m — j — 1 we 
define 

BiAt) = Bi+tj-ift) 

Mi-\-j ti 

where, as above, possible division by zero is resolved by the 
convention that a division by zero is zero. 


Theorem II.3. Let n. d be positive integers, with n > d.+ 1. Let 
Pi,... ,p n £ K 2 be n control points, and let r = 
be a nondecreasing sequence of real numbers, that is, t\ < 
... < tnfld+i- The (parametric) spline curve f(t) of degree d 


with control points {pi)™ =1 and knots r = can be 

written as: 

n 

fit ) = 'Y^PiB.flit) 

i—1 

where Bi flit) is the ith B-spline of degree d with knots r. 

The notion of linear combination of B-splines, also called 
spline function, is formalized in the following definition. 

Definition II.4. Let n. d be positive integers, with n > d + 
1. Let r = (L)”jj d+1 be a nondecreasing sequence of reals, 
and let Bi flit),B n flit) be the n B-splines of degree d 
with knots r. The linear space of all linear combinations of 
Biflit), ■ ■ ■, B n flit) is called the spline space of degree d with 
knots t and denoted by: 

^ d,r — \ ^ ^ CiBi fl | a £ M, i — 1, . . . , 71 

U=i 

An element / = YA=i c i B i,d of §d,r is called a spline function 
of degree d with knots r, and (cj)” =1 are called the B-spline 
coefficients of /. 

Now we consider the approximation problem. The method 
is based on the standard least squares approach, in which 
classically the function to minimize is the sum of squared 
errors committed at the given data points. The problem can 
be formulated as follows. 


Problem II.5. (Least Squares Approximation Problem) Let 

to > n\ given data ixu, 2/fc)fcLi an d a spline space S d, T whose 
knot vector r = 1 satisfies L+d+i > U for i = 

1 find a spline function / = X]"=i c i B iflit) £ S d, T 

which solves the minimization problem 


^iVk-fixk )) 2 = \ y k-J 2 c *Biflixi)) = 


k=l 


k= 1 


= m, in y2idk - 9ixk)f 


It is easy to see that this is a LS problem which can be 
expressed in matrix form as follows: find c = (ci,..., c rl ) f 
s.t. 



\\y- Bc \\l = 

min 

ot= (q:i ,... ,a n ) 

where 

1 

( Biflixi) 

l) 

H 

\ B 1 flAm) 

• • • -^n,d(^'ra) 


II V ~ Ba | 


e Mat„ 


( 1 ) 


is the coefficient matrix, also called B-spline collocation 
matrix, and 


Cl 


V i 


c = 


£ Mat. 


; nx i(lR), y = 


£ Mat 


•rax 11 


It is well-known that solving G is equivalent to solve the 
normal system: 

B*Bc = B f y (2) 
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Note that Problem |II.5| has always a solution. Moreover, 
there are explicit conditions under which the matrix B l B is 
invertible, and consequently system (|2|, as well as Problem 
II.5 has exactly one solution (see for instance Theorem 


5.23). 

Finally, we address the following problem: suppose we are 
given a set X = {p\ of m points in M 2 , and we 
want to construct a (parametric) spline curve C of degree d 
that approximates the points. From Problem |II.5| recalling 
Definition II.4 in order to define a spline space S d, T it is 
clear that a knot vector r is required. Obviously, different 
choices of the knot vector lead to different spline spaces and, 
consequently, to different approximating spline curves. In the 
literature (see for instance ED), there are various competing 
parametrization methods, such as: 

1) Uniform: t\ = 0 and tj = for i = 2,..., m. 

2) Chord length: t\ = 0 and tj = 
for i = 2,..., m. 

3) Centripetal: t\ = 0 and fj = U_i + 
for i = 2,..., ?7i. 


II Pi-Pi 

X?=i \\Pi+i-Pj\ 


\\pi-pi-i 


|!/2 


X?=l l|Pj + l-Pj 


, 1/2 , 


In this paper we mainly focused on the chord length method, 
though the procedure for approximating plane contours using 
spline curves (described in Section |Hl|i has been implemented 
with all three parametrization methods. 


III. Algorithm Description and Implementation 


In this section we describe an algorithmic procedure for 
approximating a given plane contour using spline curves. The 
procedure can be divided into three parts: 

I. acquisition of a digital image reproducing a silhouette; 
extraction of the contour C by strengthened edge detec¬ 
tion techniques; 

II. computation of the spline curve, approximation of the 
contour C; 

III. plot of the computed spline curve (and eventually, com¬ 
parison with C). 


Accordingly, the implementation of the algorithm is divided 
into three steps: the first and the third step, corresponding 
to parts I and III, are implemented in MatLab, using well- 
established built-in functions; the second step, corresponding 
to part II, is implemented in C++ language, using routines both 
of CoCoALib El, a GPL C++ library - the mathematical 
kernel for the computer algebra system CoCoA-5, and of the 
numerical library GSL - GNU Scientific Library fl3l . Note 
that part II is the core of the algorithm: based on Section|IIJ and 
in particular on Problem II.5 and the subsequent discussion, 
its implementation requires as input: 


11. an affine set of points in R 2 (resulting from the edge 
detection on the input image); 

12. a parameter n, corresponding to the number of polyno¬ 
mial patches; 

13. a parameter d, corresponding to the degree of polynomi¬ 
als; 

14. a parametrization method (see Section |ll]). 


It returns as output; 


Ol. a coefficient vector defining the piecewise polynomial 
function made up of n patches of degree d; 

02. the corresponding Least Square Error. 

Just to give an idea, we consider the silhouette of a horse, 
represented in Figure 1(a) Using classical edge detection 


techniques, the contour of the horse is extracted (see Figure 
1(b) I and stored as a set of 2713 points in ]R 2 . The chord 


length parametrization method is chosen (see Section [Ilj i, and 
contour approximations are obtained through computations 
with different values of the parameters, namely n = 100, d = 
2 (see Figure l(c)| >, n = 200, d = 2 (see Figure 1 1(d) i, 
77 = 100, d = 3 (see Figure 1(e) i and n = 200, d = 3 (see 
Figure We observe that there is a substantial difference in 
the approximation quality comparing the cases n = 100 and 
?z = 200, whereas we notice a slight improvement passing 
from the cases with d = 2 to d = 3. 





Fig. 1. Silhouette of a horse 


IV. Conclusion 

We have presented a novel procedure to approximate a plane 
contour by piecewise polynomial functions, depending on 
various parameters (degree, number of local patches, selection 
of knots). In order to optimize our choice of parameters, 
we plan to perform a series of behavioural experiments by 
adopting the following procedure. First, a set of contours will 
be selected, representing both shapes of daily life objects and 
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meaningless curves. Next, every contour will be approximated 
by different splines, obtained by varying parameters in such 
a way that each approximation is perceptually closed to the 
original shape. Then, participants will be asked to perform 
a same/different task. Finally, the experimental results shall 
allow us to individuate a lower dimensional spline space 
approximating the space of perceptually discriminable shapes 
in a neighborhood of a fixed shape. At this point, it will be 
possible to ask questions about shape representations such as: 

1) Is the dimensionality of the space of perceptually dis¬ 
criminable shapes the same for neighborhoods of differ¬ 
ent shapes? 

2) In particular, does the dimensionality of the space of 
perceptually discriminable shapes change as a function 
of experience? 

3) As a consequence, do neighborhoods of real objects have 
higher dimensionality than neighborhoods of meaning¬ 
less curves? 
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